Directional recoil rates for WIMP direct detection 
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New techniques for the laboratory direct detection of dark matter weakly inter- 
acting massive particles (WIMPs) are sensitive to the recoil direction of the struck 
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' nuclei. We compute and compare the directional recoil rates dR/dcos9 (where 

is the angle measured from a reference direction in the sky) for several WIMP ve- 



locity distributions including the standard dark halo and anisotropic models such 
as Sikivie's late-infall halo model and logarithmic-ellipsoidal models. Since some 
detectors may be unable to distinguish the beginning of the recoil track from its 



end (lack of head-tail discrimination), we introduce a "folded" directional recoil rate 

& . 

dR/d\cos9\, where |cos#| does not distinguish the head from the tail of the track. 

d 

We compute the CS2 and CF4 exposures required to distinguish a signal from an 
■ isotropic background noise, and find that dR/d\cos9\ is effective for the standard 

m . 

, dark halo and some but not all anisotropic models. 

o 
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CSj . PACS numbers: 95.35.+d 
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I. INTRODUCTION 

The nature of dark matter (DM) in the Universe is still one of the outstanding problems in 
astrophysics and cosmology. Numerous observations support the existence of DM. Examples 
are: big bang nucleosynthesis, cosmic microwave background data (WMAP3), supernova 
surveys, Galaxy surveys (SDSS, 2dF), and distance measurements with Cepheids (HST). 
In the concordance cosmological model ACDM, the total density of the Universe has three 
contributions: matter, radiation, and a cosmological constant. The matter contribution can 
be further divided into the contribution of ordinary (baryonic) matter and the contribution of 
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non-baryonic cold dark matter (CDM). Their density parameters, i.e. their densities in units 
of the critical density p c = 1.053 x 10~ 5 /i 2 (GeV/c 2 ) cm~ 3 (here h is the Hubble constant in 
units of 100 km/s/Mpc), are: fi b = (0.02186 ± 0.00068) h' 2 and Q CDM = (0.1105±g;gg||)^- 2 
(from Ref. [1]). Thus, CDM constitutes ~ 84% of the matter in the Universe. 

CDM is found in clusters of galaxies and in individual spiral and elliptical galaxies. For 
example, stars in spiral galaxies are observed to move too fast around their galactic centers to 
be explained by the gravity of luminous matter alone. In particular, our Milky Way Galaxy 
also contains DM. Binney and Dehnen [2], for example, show that the rotation curve of the 
Milky Way is nearly constant far beyond the Sun's location, implying the presence of DM 
in the Sun's neighborhood. This and similar studies give a density of DM near the Sun of 
p = 0.3 (GeV/c 2 )/cm 3 . 

The nature of CDM is still unknown. DM candidates for CDM are subatomic particles 
such as neutralinos, axions, Kaluza-Klein particles, and other WIMPs. WIMPs are hypo- 
thetical electrically-neutral stable particles with scattering cross-section off nucleons of the 
order of the weak interaction (<r p 10~ 44 cm 2 ) and mass in the range 10 — 1000 GeV. Dark 
matter WIMPs arise for example as lightest supersymmetric particles (LSPs) in supersym- 
metric extensions of the Standard Model of particle physics. Kaluza-Klein particles arise 
in theories with more than four space-time dimensions, and share the same properties with 
WIMPs except for being somewhat heavier. In this paper, we refer to WIMPs but our 
considerations apply to Kaluza-Klein particles as well. 

Dark matter WIMPs near the Sun can reach the Earth and can scatter elastically off 
target nuclei in a detector, making the nuclei recoil. The energy, and recently the direc- 
tion, of the recoiling nuclei can be measured experimentally. Extensive experimental efforts 
have been devoted to detect WIMPs directly (e.g., DAMA, CDMS, EDELWEISS, CRESST, 
DAMA/LIBRA, SuperCDMS, DRIFT, etc.). There are two types of direct detection ex- 
periments: those that measure the recoil direction and those that do not. Examples of 
directional direct detectors are: DRIFT [3-7], which uses a gas target in a time projection 
chamber and has run a prototype detector for a few years; NEWAGE [8, 9], which uses a 
similar time projection chamber and is sensitive to spin-dependent WIMP-nucleus interac- 
tions; and detectors that use organic crystals such as stilbene [10, 11]. The other detectors 
previously listed are all non-directional. 

A goal of directional WIMP detectors is to identify Galactic WIMPs by using the dis- 
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tribution of the nuclear recoil directions as a signature. We believe that an analysis of the 
WIMP-induced recoil directions can also allow the study of the structure and dynamics of 
the WIMP halo. 

The idea of directional WIMP detection originated as early as 1988. Spergel [12] suggested 
that a WIMP signal could in principle be identified by means of the diurnal rotation of the 
'WIMP wind' direction due to the Earth's rotation (the 'WIMP wind' is caused by the Solar 
System's rapid motion through the Galactic halo). The practical realization of directional 
WIMP detection was delayed by the difficulty of finding a suitable target material and an 
effective detection technique. In 1996, Martoff et al. [13] described a prototype direction- 
sensitive solid-state detector for WIMPs. A gaseous directional detector was studied by 
Martoff et al. [3] and is described in Snowden-Ifft, Martoff, and Burwell [4] (see also [5-7]). 
This detector, called DRIFT (Directional Recoil Identification From Tracks), uses a time 
projection chamber filled with a low pressure mixture of a target gas and an electronegative 
gas. The first stage of DRIFT (DRIFT I) had a 1 m 3 target (167 g of CS 2 ) and ran from 2001 
to 2004 at the Boulby mine, North Yorkshire, England [14]. The current stage of DRIFT 
(DRIFT II) is an array of 1 m 3 modules and has been operational since 2005 [15]. A future 
stage has been envisaged (DRIFT III) that may have a target mass of up to 1 ton [16]. 

On the theoretical side, Copi, Heo, and Krauss [17] and Copi and Krauss [18] examined 
the number of events required to distinguish a WIMP signal from an isotropic background. 
Gondolo [23] obtained analytic expressions for a variety of directional recoil spectra by means 
of the Radon transform that relates the WIMP velocity distribution to the distribution of 
recoil momenta. Freese, Gondolo, and Newberg [25] studied the possible directional detection 
of WIMPs belonging to the Sagittarius tidal stream, which may be showering DM onto the 
Solar System. Morgan, Green, and Spooner [26], Morgan and Green [27], and Green and 
Morgan [28] studied how the exposure required to directly detect a WIMP directional recoil 
signal depends on the capabilities of a directional detector. They also examined statistical 
tests to distinguish a WIMP signal from an isotropic background and found that in detectors 
with head-tail discrimination (see below) of order ten events will be sufficient to distinguish 
a WIMP signal from an isotropic background for all of the halo models they considered. 
Host and Hansen [29] investigated the possibility of measuring the velocity anisotropy of 
the Galactic dark matter halo in a direction-sensitive WIMP detector. They found that 
in excess of 10 5 events across all energies are needed to make a coarse measurement of the 
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velocity anisotropy. 

The goal of this paper is to study how different halo models affect the directional recoil 
rate dR/dcos9, where 9 is the angle between the nucleus recoil direction and a chosen 
reference direction in the sky. We compare the directional recoil rates for the different 
models. We repeat the same analysis for a "folded" directional recoil rate dR/d\cos9\ that 
incorporates the inability of some detectors to distinguish the beginning of a recoil track from 
its end (head-tail discrimination). We compare each dR/d\cos9\ to an isotropic background, 
to examine the possibility of discriminating a WIMP signal from background noise. 

In Sec. II, we present a general discussion of the directional recoil spectra. There we give 
the expressions of various differential recoil rates that are useful to analyze and interpret 
WIMP direct detection experiments. In Sec. Ill, we describe two methods, numerical and 
analytical, for calculating the directional differential recoil rate dR/dcos9 of recoiling target 
nuclei struck by WIMPs. The analytical method is applied to a Gaussian velocity distribu- 
tion whose average velocity is aligned with the reference direction. The numerical method 
is more general and it can be used for any reference direction and any WIMP distribution. 
In Sec. IV, we present the results of applying the numerical method to various WIMP halo 
models, including streams of WIMPs, the standard dark halo, the Sikivie late-infall halo, 
and anisotropic models. In Sec. V, we address the difficulty of head-tail discrimination 
in WIMP direct detection experiments and present recoil distributions suitable for direct 
comparison with experiments lacking head-tail discrimination. Finally, we summarize our 
results in Sec. VI. 

II. DIRECTIONAL RECOIL SPECTRA 

In this section we give an expression for the directional recoil rate for interactions between 
WIMPs and target nuclei. In WIMP direct detection, the collision between the WIMP and 
the target nucleus is detected by measuring the energy of the recoiling nucleus. In directional 
detectors, one can also measure its direction of recoil. 

Fig. 1 shows the kinematics of such a collision. The energy of the recoiling nucleus is 
given by (see e.g. Gascon [30]) 

E = £ max cos 2 R , (f) 
where 9 R is the angle of the nuclear recoil relative to the initial WIMP direction (recoil 
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FIG. 1: Kinematics of a WIMP-nucleus elastic scattering. 

angle), and 

F - (9\ 

max " M n ' [) 

is the maximum energy that the WIMP can transfer to the nucleus. Here v is the speed 
of the incoming WIMP, m is its mass, M n is the mass of the target nucleus, and /j, n = 
m M n j (m + M n ) is the reduced mass of the WIMP-nucleus system. 

In general, the differential recoil spectrum, i.e. the differential event rate per unit detector 
mass, is given by 

where the sum is over the nuclear species in the target, C n is the fraction of mass in species 



is the minimum WIMP speed required to transfer an amount of energy E to the nucleus of 
mass M n in the detector (here c is the speed of light), p is the local WIMP density mentioned 
in the Introduction, £{E) is the detection efficiency at recoil energy E, and a n (E) is defined 
as 

a n (E) = E mK ^ (5) 
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FIG. 2: The figure shows the reference direction n, pointing to a specific direction in the sky, and 
the angles 6 and (f> used in Sec. III. 

with da/dE equal to the differential WIMP-nucleus scattering cross section. 

For directional detectors, we need a differential rate not only in energy but also in direc- 
tion. The three-dimensional recoil rate in spherical coordinates where the angles 9 and 
refer to the direction of the nuclear recoil and the radial coordinate is the recoil energy E, 
is given by (see Ref. [23]) 



is defined to be the integral of /(v) on a plane orthogonal to the direction w at a distance 
w from the origin v = [24]. In formulas, 



where 5 is the Dirac delta function and w is the recoil direction. In this work we will specify 
the direction w using a reference frame fixed in the sky in preference to a reference frame 
fixed with the laboratory. 





where dQ = d(j)d cos 9 and f(w,w) is the 3-dimensional Radon transform of the velocity 
distribution function /(v). The 3-dimensional Radon transform f(w,w) of a function /(v) 




(7) 
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Projections of the directional differential rate are also useful and have been used in the 
past. For example, one can measure recoil directions w from a chosen reference direction 
n as in Fig. 2. If Eq. (6) is integrated over the azimuthal angle <p and the energy E, one 
obtains 

dR f f dR ..^ (g) 



dcos9 J J dEdVt 

where 9 is the angle between the reference direction n and the recoil direction w. Eq. (8) 
is the directional differential recoil rate we study in this paper. It has been used in previ- 
ous work [8, 31] to compare WIMP velocity distributions and/or assess the advantages of 
directional detection methods. 

The directional rate dR/dcos9 requires a 3D read-out of the track direction. Although we 
are optimistic that one day a 3D read-out will be available, current experiments are limited 
to a 2D read-out in a plane fixed with the laboratory [15]. This plane precesses around 
the North-South terrestrial axis due to the rotation of the Earth. A differential rate dR/d<p 
appropriate for this situation has been introduced and studied in Refs. [26-28], to which we 
refer. 

Besides the difficulty of a 3D read-out, current detectors may be unable to distinguish the 
beginning of the recoil track (the head) from the end of the track (the tail). This is called the 
difficulty of head-tail discrimination. Because of this, it is useful to introduce the following 
"folded" directional recoil rate relevant to experiments that lack head-tail discrimination: 

dR _dR(cos9) dR(-cos9) 

H TZIZn • ( 9 ) 



d\cos9\ dcos9 dcos9 
This rate is correctly normalized because the integral of both sides gives the total rate. As 
illustrated in Fig. 3, |cos#| does not distinguish between w and — w. Therefore, there is no 
need to know the track's head from its tail when using the "folded" directional recoil rate, 
Eq. (9). However, there may be a loss of information in doing so (see Section V). 

The differential WIMP-nucleus scattering cross section a n (E) in Eq. (5) can be split into 
two parts, one spin-independent (SI) and the other spin-dependent (SD): 

a n {E) = a s n I {E) + a s n D {E). (10) 

Correspondingly, one can separate the spin-independent and spin-dependent contributions 
to the directional recoil rate dR/d cos 9 as 

'/ 1 { (i n 1 (i n j (li) 



d cos 9 d cos 9 d cos 9 
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FIG. 3: The vectors w and — w share the same value of |cos#|. 

A similar separation can be defined for the "folded" directional recoil rate dR/d\cos9\. 

The rest of this section describes the expressions for the directional recoil rate for spin- 
independent and spin-dependent interactions. 

A. Spin-independent Directional Recoil Rates 

In Eq. (10), the spin-independent part a^ T (E) can be written as 

a s n I (E)=a F n (E), (12) 

where <7 is the WIMP-nucleus scattering cross section and J r n (E) is a nuclear form factor 
which depends on the type of WIMP-nucleus interaction and on the mass and spin distribu- 
tions within the nucleus. In cases where the nuclear form factor effects are negligible we have 
J- n {E) — 1. In reality, the nuclear form factor may become important for specific detectors. 
One can write 

a = ^\ZG p s + {A-Z)G n s \ 2 , (13) 

7T 

where Z is the number of protons in the nucleus, A is the mass number of the nucleus, and 
(jP (GJ) is the effective proton (neutron)-WIMP coupling. The WIMP-proton cross-section 
is 

^p = ^W, (14) 

71 
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where 



Hp = : (15) 



m + m p 



is the WIMP-proton reduced mass. Assuming, as it is approximately the case for neutralino 
dark matter, that 

g p s = g:, (i6) 



we have 



a = ^A 2 a p . (17) 



In this case, the recoil rate, Eq. (6), takes the form 
dR SI po p 



f n 

We can define an effective spin-independent recoil momentum distribution, f^(E,w), as 
the average over all masses 

fS(E, w) = J2Cn A 2 n f(w n , w) F n {E) £{E) . (19) 



The rate of detection of WIMPs then reads 

dR SI pa p z SI 



dE dfl 4 7r fi p m 

We can also write 



f%(E,*). (20) 



dcos 9 4 7i /Xp m 
For example, for a CS 2 target as in DRIFT, 



^ 1 1 fg(E,*)d4>dE. (21) 




° S ~ 2M S + M C ' (22) 
° C = 2M^+M C ' (23) 

and fcff(E, w) is given explicitly by 

/cff(E ' W) " 2 .U s • Mr £ - (24) 

Here f n = f(w n ,w). (Notice that the symbol C (italic) denotes the fraction of mass while 
the symbol/subscript symbol C (roman) denotes the carbon nucleus). 
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Using common units and magnitudes, the spin-independent directional detection rate of 
WIMPs is 

£bn = 1-306 x 10- 3 - fcs(E, w) , (25) 

dEdVt kg-day-keV-sr ^ \i\m JeSK ' n K J 

where p .3 is the DM density in the solar neighborhood in units of 0.3 (GeV/c 2 )/cm 3 , 044 

is the proton cross-section in units of 10~ 44 cm 2 . fx p and m are in GeV/c 2 , and /J/ is in 

(km/s)- 1 . 

B. Spin-dependent Directional Recoil Rates 

In Eq. (10), the spin-dependent part a^ D (E) can be written as 

QO ,.2 fl2 

a * D{E) = ( 2 i + # [ °p + a ° + ap ° n • (26) 

Here h is the reduced Planck constant, Gf is Fermi coupling constant (GV/(^-c) 3 = 1. 16637 x 
10~ 5 GeV~ 2 ), J n is the nucleus total angular momentum in units of H, a p (a n ) is the effective 
axial coupling of WIMP and proton (neutron) in units of 2y/2G F /h 2 [19]. (Notice that the 
subscript/superscript symbol n (italic) denotes the nucleus while the subscript symbol n 
(roman) denotes the neutron). In Eq. (26), the dimensionless functions S pp (E), S™ n (E), and 
S™(E) play the same role as the nuclear form factor T n [E) in the spin-independent case. 
They are given by 

Spp = S5, + S?i + S n m , (27) 

nn qn , qn qn /no"! 

°nn _ °00 °11 °01 J l ZC 7 

5^ = 2(55, - £&) , (29) 

where Sq , S^, Sq X are the nuclear spin structure functions defined in [20]. When the nuclear 
spin is approximated by the spin of the odd nucleon only, one finds 

S; p = A " Jn(J " + 1)(2Jn + 1) , = 0, S; n = 0, (30) 

7T 

for a proton-odd nucleus, and 

rin _ n on _ \i Jn(Jn + l)(2J ra + 1) on _ n /oi \ 

PP ~~ ' nn _ ) °pn — U ) 

71 

for a neutron-odd nucleus. Here A n is conventionally defined through the relation <n\S\n> = 
\ n <n\ J\n>, where \n> is the nuclear state, S is the spin, J is the total angular momentum. 
Tables of A 2 J n (J n + 1) values for several nuclei can be found in [21] and [22]. 
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The spin-dependent cross-section off a proton is 



24fi 2 p G 2 F 

% =—^ a ^ ( 32 ) 



and that off a neutron is 

sd _ 24/i n (jp 2 



a 2 n . (33) 



In case the target is a combination of different nuclei, we write 

§5S = S K *> + a » *> + a ' A *>] ' (34) 

where we define the effective spin-dependent recoil momentum distributions 

/Sp(£, w) = 3{2 j* + 1} Cn S; p (E) f(w n , w) 5(£'), (35) 

and similarly for /J F D nn (£', w) and / c ff D pn (-E, w). The normalization of f^ PT ,(E, w), 
/eff D nn(^>w), and /eff D pn (£,w) has been chosen so that f^ s D pp {E,w) = f{w p ,w) when the 
target is a proton and f^ njl (E, w) = /(u> n ,w) when it is a neutron. 

For the CF4 target used in the NEWAGE detector [8, 9], and in the proton-odd approx- 
imation, Eq. (34) takes the form 

In this case, the C nucleus has no spin (thus Ac = 0), while the F nucleus has spin |. In the 
proton-odd approximation, A|Jf(Jf + 1) = 0.647 (see Table 1 in [21]). Thus the effective 
spin-dependent recoil momentum distribution, Eq. (35), reads 

/Sp( CF 4) = ^ 0.647 C F f F S(E). (37) 
The fraction of mass C F is given by 

C F = 4Mf . (38) 

F 4M F + M C 1 ; 

Using common units and magnitudes, the spin-dependent directional detection rate of 
WIMPs (for the CF 4 target used in the NEWAGE detector) is 

dRSD _ 1 ana v m-3 eVelltS ^0-3 <44 fSD (p ~x ,o Q N 



dE dVt kg-day-keV-sr 4 7r fi p m 
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III. CALCULATION OF dR/dcos6 



The analysis of directional spectra in WIMP direct detection can be carried out by 
computing the directional differential recoil rate dR/dcos6 as a function of the angle 9 
between the nuclei's recoil directions w and a reference direction n. In this section we 
describe two methods for calculating the directional recoil rate dR/dcos8. The first method 
is numerical and can be used for any WIMP velocity distribution and for any reference 
direction. The second method is analytic and is restricted to Gaussian distributions and to 
reference directions n in the same direction as the WIMP average velocity V. Results of 
applying these methods to various dark halo models are given in Sec. IV. 

Before describing the two methods, we recall the expression of the recoil momentum 
function f(w, w) for a Maxwellian velocity distribution. In the rest frame of the detector, 
the WIMP velocity distribution is given by 



where v is the velocity of a WIMP, a v is the velocity dispersion (not to be confused with the 
WIMP-nucleus cross section) and V is the average velocity of the WIMPs with respect to 
the detector. The recoil momentum spectrum for nucleus n in the laboratory frame is [23] 



where V(W, G) is the average velocity of the WIMPs relative to the Galactic rest frame 
(zero in the standard halo model), V(S, G) is the velocity of the Sun relative to the Galactic 
rest frame (of order 200 km/s), V(E, S) is the velocity of the Earth relative to the Sun (of 
order 30 km/s), and V(lab, E) is the velocity of the detector in the laboratory relative to 
the center of mass of the Earth (of order 0.3 km/s). As the Earth rotates and orbits the 
Sun, two signal modulations (annual and diurnal) are expected as a result of the relative 
motions V(E, S) and V(lab, E). In this work we neglect V(E, S) and V(lab, E) - and the 
corresponding annual and diurnal modulations - and use the velocity distribution in the 
frame of the Sun. We also neglect focusing effects due to the gravitational field of the Sun 
(of order 1 km/s, see [32]). 




(40) 




(41) 



In principle, 



V = V(W, G) - V(S, G) - V(E, S) - V(lab, E), 



(42) 
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In the following, we assume that the nuclear form factors are T n [E) — 1. In the case of a 
detector with threshold energy E thr , we model the detection efficiency at recoil energy E as 



In general directions and for non-Gaussian distributions, the integration of the recoil 
rate over the energy E and the azimuthal angle in Eq. (8) cannot be done analytically. 
However, one can calculate the recoil distribution numerically. 

We do the integration over the azimuthal angle <fi using a Riemann sum. For the inte- 
gration over the energy E, it is difficult to use a Riemann sum because of narrow peaks 
in dR/dcos9 as a function of cosO when streams of dark matter are present. Therefore, 
we do the integration over E by means of the fifth-order Cash-Karp Runge-Kutta method 
with adaptive stepsize control, as described in Ref. [33]. To apply this method, we write the 
differential recoil rate in the form of an ordinary differential equation 



where A0 = 2n/N ( f > , and the unit vectors w i; which depend on fa, are specified in Eq. (52) 
below. We take = 100 with initial condition (dR/d cos 0) e=o = 0; we require an accuracy 
of 10~ 8 and choose a scaling factor 



where /irk is the value of the stepsize, yi = dR/dcosO at = fa, and (dy/dE)i is the partial 
derivative of dR/dcos9 with respect to E at = fa. 

An expression for Wj is obtained as follows. At a fixed value of the angle between 
the reference direction n and the recoil direction w, the possible directions of the target 
nucleus recoils lie on a cone of opening angle 20, as shown in Fig. 2. Since is constant, it 
is convenient to specify w in terms of the polar angle and the azimuthal angle reckoned 
with respect to the n axis. 

We want to find the Cartesian components of Wj in an arbitrarily given reference frame 
x, y, z. For this purpose, we introduce an auxiliary Cartesian coordinate system x', y', z' , 




(43) 



A. Numerical dR/dcos - general case 




(44) 




(45) 
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with z' aligned with n. Since 9 and are reckoned from z' , the vector Wj, which is w at 
4> = fa, can be written in terms of the new unit vectors e^, e' 2 , e' 3 as: 

Wj = sin 9 cos fa e[ + sin 9 sin 0j e' 2 + cos e' 3 . (46) 

Our task is now to express the new basis vectors e^, e' 2 , e' 3 in terms of the original basis 
vectors e?i, e 2 , e 3 . Once this is achieved, Eq. (46) will give Wj in terms of ei, e 2 , e 3 and one 
can read off its Cartesian components. 

The transformation matrix from ei, e 2 , e 3 to e^, e' 2 , e' 3 can be found using Euler angles. 
We first rotate the ej axes about the z axis counterclockwise through an angle a + |. Then 
we rotate the resulting axes about the new x axis counterclockwise through an angle (3. We 
find 

= — sin a ei + cos a e 2 , (47) 
e' 2 = — cos j3 cos a §i — cos (3 sin a e 2 + sin (5 e 3 , (48) 
e' 3 = sin P cos a e x + sin (3 sin a e 2 + cos (3 e 3 . (49) 

Identifying n with e 3 in Eqs. (47-49) shows that a and (3 are the spherical coordinates of 
n: 

a = tan' 1 (^j , (50) 
(3 = cos' 1 (n z ) . (51) 

Inserting these relations into Eq. (46), we obtain 

Wj = (— sin 9 cos 0« sin a — cos (3 cos a sin 9 sin fa + sin (3 cos a cos #) §i 
+ (sin ^ cos 0j cos a — sin 9 sin 0j cos /3 sin a + sin /3 sin a cos ^) e 2 
+ (sin 9 sin 0j sin (3 + cos /3 cos 9) e 3 . (52) 

B. Analytic dR/dcos6 - special case 

We were able to derive an analytic expression for dR/dcos9 when the WIMP velocity 
distribution is Maxwellian and the reference direction n is aligned with the WIMP average 
velocity V. Notice that this does not mean that 9 = 9 R , the recoil angle (see Fig. 1), 
because 9r is measured from the velocity of an individual WIMP while 9 is measured from 
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the average velocity of all WIMPs. Here we assume a spin-independent case and a zero 
threshold detector. 
We wish to compute 

dR pa p 



Using Eq. (19), 



where 



d cos 9 4 7r /J,pTn 
dR 




d cos 9 4 7r fx 2 m 



f eS (E,w) dEd<p. 



^1 CnA^Im 



poo plix 

I n = / f n (w,w)dEd(f). 
Jo Jo 



For a Maxellian velocity distribution, 

■y poo p2n 



K-w-v] 2 ', ,,. , 

' S 1 } 2a 2 )dEd(j). 



(53) 



(54) 



(55) 



(56) 



Since w ■ n = cos 9 and we assume here that V = Vh, we have 



w • V = V cos 9 . 



(57) 



From Eq. (4), 



4/i 2 

dE = n 9 w„. dw r . 
M n c 2 



Substituting Eqs. (57) and (58) into Eq. (56), we find 

87T/4 f°° _ / [w n -Vcos9] 2 



n ' M n c 2 (2na 2 y/ 2 



cxp 



2 at 



M n c 2 (2tt)V2 



cxp 



V 2 cos #\ Ucosi 



2 a 2 



w n dw n , 



1 + erf 



Vcos9 



a v V2 



where erf(x) is the error function. Inserting Eq. (60) into Eq. (54) gives 



dR 2p<J v C n A 2 n n 2 n o v ^ \/ 2 cos 2 ^ + Vcosfl ^- 



dcos# ^ V2nM n c 2 mfil 



2al 



a v \/2 



. Vcos9\ 
1 + erf ( pr- 



or 



rf-R _ 2612 events po.3cr 44 C n A 2 n iJ, 2 n a v 
dcos9 kg-day ^ v / 27rM n c 2 m/i 2 

V 2 cos 2 #\ Vcostf 
exp | — ) H —\/tt 



2a 2 



a v V2 



1 + erf 



/Vcos9 

V (T v y/2 



(58) 

(59) 
(60) 



(61) 



(62) 



This analytic formula agrees with the numerical method, Eq. (44), in the cases in common. 
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IV. RESULTS 

In this section we present the results of applying the numerical method, Eq. (44), to 
various WIMP halo models, including streams of WIMPs, the standard dark halo, the 
Sikivie late-infall halo (SLI streams), and anisotropic logarithmic-ellipsoidal models. For 
comparison's sake, in the last three cases (the standard dark halo, SLI streams, and the 
anisotropic models), we fix the reference direction n to be the direction of Galactic rotation. 
For all these models, we use CS 2 molecules as target nuclei. We use ecliptic coordinates 
(longitude A and latitude (3). 

We specify the WIMP velocities and the directions of nuclear recoil in a Cartesian coordi- 
nate system defined as follows: as seen from the Earth, the x axis points toward the position 
of the Sun at the vernal equinox, the y axis toward the position of the Sun at the summer 
solstice, and the z axis toward the North Pole of the ecliptic (which is the projection of 
the trajectory of the Sun onto the celestial sphere). The ecliptic longitude of the reference 
direction, A ax i S , is the angular distance along the ecliptic from the vernal equinox to the 
base of the great circle containing n and the pole of the ecliptic; it is measured eastwards in 
degrees from 0° to 360°. The ecliptic latitude of the reference direction, /3 axis , is the angular 
distance north (from 0° to 90°) or south (from 0° to —90°) of the ecliptic along the previously 
mentioned great circle; it is measured from the ecliptic to n. In terms of A axis , /3 ax i S , we write 

n= (cos/3 axis cosA axis , cos/9 axis sinA axis , sin/3 axis ) . (63) 

In computing our rates we take a WIMP mass m = 60 GeV, a WIMP-proton cross section 
a p = 10~ 44 cm 2 = 1CT 8 pb. We also neglect the nuclear form factor, i.e. take F n (E) = 1. 
In this section, we further assume an ideal efficiency with zero threshold S(E) = 1. In 
Sec. V, we present some results with a non-zero threshold for both spin-independent and 
spin-dependent cases. 

A. Streams of WIMPs 

We start by showing the results of applying the numerical method, Eq. (44), to a simple 
model, namely a Maxwellian stream of WIMPs with an average velocity V. For definiteness 
sake, we consider the case of a stream with V = (V, 0, 0), but our results apply to a generic 
V. 
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FIG. 4: The directional recoil rate dR/dcosO off a CS2 target as a function of cos# for streams 
with average velocity V parallel to the reference direction n for different a v /V ratios; a v /V = 
4.74/300 = 0.0158 (solid line), a v /V = 10/100 = 0.1 (dashed line), a v /V = 100/100 = 1 (dotted 
line), and a v /V = 200/100 = 2 (dashed-dotted line). 

First we fix the reference direction n in the direction of V, i.e. we take (Aaxis, /3 ax is ) = 
(0°,0°), and vary the ratio <r v /V of the velocity dispersion a v to the magnitude of the 
average WIMP velocity V. Fig. 4 shows the resulting dR/dcos9 as a function of cos 6* for 
streams with a v /V = 4.74/300 = 0.0158 (solid line), a v /V = 10/100 = 0.1 (dashed line), 
a v /V = 100/100 = 1 (dotted line), and a v /V = 200/100 = 2 (dashed-dotted line). 

The solid line in Fig. 4 shows that dR/dcos8 peaks at 9 = and almost vanishes for 
cos# < 0. This behavior can be understood by considering a stream with zero velocity dis- 
persion. All WIMPs in the stream move at the same velocity V, so the velocity distribution 
function is 

/(v) = *(v - V) . (64) 
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Its Radon transform follows from Eq.(7) as 

f(w,w) = J S(v w -w) f(v)d 3 v = S(V w -w), (65) 
With V = Vh and w • n = cos 9, 

f(w,w) = 5(w-Vcos9). (66) 
Inserting this relation into Eq. (55), and using Eq. (58), gives 

In = j I ^ w ^ dEd< ^ ( 67 ) 

poo 

= 2ix j 5(w n - V cos 9) dE (68) 
Jo 

POO 

2n / S(w n — V cos 9)w n dw n (69) 
Jo 



M n c* 

8 S vms<> - (70) 



Therefore, for zero velocity dispersion, we have 



dcos9 




2 pa v> ji 2 n C n A 2 n 

filmM n c 2 J ' (71) 
0, ifcos^<0^ 



For positive cos#, dR/dcos9 is a linearly increasing function of cos#. Its maximum occurs 
in the forward direction at 9 = 0°. Away from the forward direction, the number of recoils 
decreases, and dR/dcos9 drops. At negative cos# there are no recoils for a v = 0, since 
momentum conservation forces all recoils to be in the forward direction. 

As the ratio a v /V increases we can observe first a few and then many recoils at cos 9 < 0, 
because of the effect of the relatively higher velocity dispersion a v of the streams. This is 
illustrated by the dashed, dotted, and dashed-dotted lines in Fig. 4, some of which extend 
to cos# — —1. 

Now we fix the velocity of the stream and its dispersion, and vary the reference direction 
n. Since there are no other directions in this case, the rate dR/dcos9 depends only on the 
angle ip between V and n, besides the ratio a v /V. For definiteness, we take V = (200, 0, 0) 
km/s and a v = 10 km/s (so the ratio a v /V = 0.05). We start with the stream's average 
velocity V parallel to the reference direction (A axis , /9 axis ) = (0°,0°). We keep the ecliptic 
latitude of the reference direction constant at /3 axis = 0°. We increase its ecliptic longitude 
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FIG. 5: The directional recoil rate dR/dcos9 off a CS2 target for streams with V = 200 km/s, 
a v = lOkm/s (thus a v /V = 0.05) and various reference directions making an angle ip with V equal 
to: 0° (solid line), 45° (dashed line), 90° (dotted line), 135° (dashed-dotted line), and 180° (thick 
solid line). 

A axis in steps of 45° until V is anti-parallel to n. Fig. 5 shows the corresponding dR/dcos9 
for i) = 0°,45°,90°, 135°, and 180°. We note the following. 

First, the directional recoil rate dR/dcos6 of streams parallel to the reference direction 
peaks in the direction opposite to the incoming WIMPs. For example, the solid line that 
peaks on the right is for WIMPs coming from the left and the thick solid line that peaks on 
the left is for WIMPs coming from the right. 

Second, for a given stream, if we take two reference directions that are opposite to each 
other, they form angles ip and n — ip, respectively, with the stream's velocity V. Since 
cos^ = — cos(7r — ■?/>), their respective dR/dcos6 distributions transform into each other 
under the substitution cos 6 — > — cos 9. For example, the case -0 = 0° (solid line in Fig. 5) is 
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the reflection about cos 9 = of the case tp = 180° (thick solid line in Fig. 5). Similarly, the 
case ip = 45° (dashed line in Fig. 5) is the reflection of the case ip = 135° (dashed-dotted line 
in Fig. 5). For the same token, the case ip = 90° is the reflection of itself, i.e. it is symmetric 
under cos 9 — > — cos 9. 

B. Standard dark halo 

Here we consider a flow of WIMPs according to the standard dark halo. In this model, 
the WIMPs are on average at rest relative to the Galaxy, and their velocity distribution is 
Maxwellian with velocity dispersion given by 

220 km/s 

d = — 71 — ' 

The Local Standard of Rest (LSR) moves at t> L sR — 220 km/s relative to the Galactic rest 
frame in the direction of the Galactic rotation, i.e. (/cai.rot.j ^Gai.rot.) = (90°, 0°) in Galactic 
coordinates and (A Ga i.rot. ? Asai.rot.) = (347?340, 59?574) in ecliptic coordinates. 

Using the ecliptic Cartesian coordinate axes defined in the beginning of the current sec- 
tion, the average WIMP velocity components with respect to the Sun are 

V std = (-104.525, 34.947, -196.836) km/s . (73) 

The average velocity V st d points in the direction (A st d, Pstd) = (161?513, — 60?755). These 
values (V st( j and its direction) are from Alenazi and Gondolo [32]. 

In this subsection, we choose the reference direction n along V std , i.e. (A axis , /3 axis ) = 
(161?513,-60?755). 

Applying the numerical method, Eq. (44), we see in Fig. 6 that dR/dcosO is an increasing 
function of cos#. It peaks in the forward direction (9 = 0°) because most of the recoils occur 
at 9 = 0°. Away from the forward direction, dR/dcos9 decreases because fewer recoils result 
due to momentum conservation. The effect is similar to the dotted line in Fig. 4, whose 
ratio <y v /V = 1 is close to the ratio cr vs td/Vstd = = 0.689 of the standard dark halo. 

We further discuss the case of the standard dark halo (Fig. 6) in the following subsections. 
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FIG. 6: The directional recoil rate dR/dcosO off a CS2 target for the standard dark halo for a 
reference direction n in the direction (A ax i s , Aixis) = (A s td,Astd) = (161?513, — 60?755). 

C. Sikivie's late-infall halo model (SLI streams) 

Here we consider Sikivie's late-infall (SLI) halo model [34-36]. The SLI model is a self- 
similar axially symmetric infall model with net angular momentum and parameters adjusted 
to describe our Galaxy. In this model, many flows of collisionless DM particles are oscillating 
into and out of the Galaxy in pairs. The first pair corresponds to particles passing through 
the Solar System from opposite sides of the Galactic Plane for the first time, the second pair 
corresponds to particles passing for the second time, and so on (from the fifth pair onward, 
the flows in a pair come on the Galactic Plane but one inwards and the other outwards). 
Table I lists the densities pi and velocities Vj of the first 20 pairs of SLI streams in our 
ecliptic coordinate system. These values are taken from Ref. [37], where they were given in 
the Galactic coordinate system (see also [38]). 

The conversion from the Galactic frame to ecliptic frame proceeds as follows. In the 
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TABLE I: Densities and velocities of the first 20 pairs (±) of SLI streams in the ecliptic coordinate 
system. The values are from Ref. [37] where they were given in the Galactic coordinate system. 
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CIO -26 ff/cm 3 ") 


( 1cm / si 


( 1cm / ^ 1 


( 1cm / ^ 1 


1± 
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-560 
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-377 
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1.0 


-417 
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6.8 
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-270 
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331 


85 


8± 


1.4 
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1.1 
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-30 


-339 


-67 


8.8 
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-337 


-84 


-1.0 


369 


-16 
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0.8 


-50 
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-34 
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0.7 


-57 


-323 


-114 


-19 


363 


-47 
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0.6 


-64 


-316 
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-27 


359 


-61 


15± 
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-68 
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-70 
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0.55 
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-299 
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-37 


347 


-79 
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0.5 


-78 


-293 


-150 


-42 


343 


-88 


18± 


0.5 


-80 


-283 


-153 


-45 


334 


-94 


19± 


0.45 


-82 


-277 


-157 


-48 


329 


-98 


20 ± 


0.45 


-84 


-271 


-161 


-51 


325 


-103 



Galactic coordinate system, X points toward the Galactic Center, Y toward the direction 
of Galactic rotation, and Z toward the North Galactic Pole. Our ecliptic coordinate system 
assumes x pointing toward the vernal equinox, y toward the summer solstice, and z toward 
the North Pole of the ecliptic. We take the Galactic components of the solar motion to be 
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[39] U = 10.00 ± 0.36 km/s (radially inwards), V = 5.25 ± 0.62 km/s (in the direction of 
Galactic rotation) and W = 7.17±0.38 km/s (vertically upwards). The velocities Vf of the 
SLI streams relative to the Sun follow, in Galactic coordinate system, as 

Vf x =Vg-U, (74) 

V& = Vg-V- v LSR , (75) 

V& = Vg-W, (76) 

where Vf are the velocities relative to the Galaxy (and t> L sR — 220 km/s, as in the previous 
subsection). Then we convert Vf x , , and Vf z from the Galactic coordinate system to our 
ecliptic coordinate system, and obtain the velocities of SLI streams V ix , V iy , and V iz listed in 
Table I. 

In the SLI model, the WIMP velocity distribution function is given by 

/sli(v) = - Pi^v - V,) . (77) 

P i 

Following [38], we assume a Gaussian distribution of velocities for each flow with velocity 
dispersion <7„sli = ^ with vq = 30 km/s. We have 

/su(v) 4?(^£^ exp (-^f)' (78) 

Now, the recoil momentum distribution is 

f I M 1 Pi ( [^-W-Vi] 2 \ 

^'"'^P^yr 2.; SL1 J- < 79 » 

We fix the reference direction opposite to the direction of Galactic rotation, (/axi S , 6 ax i S ) = 
(270°, 0°) in Galactic coordinates or (A^, /3 axis ) = (167?340, — 59?574) in ecliptic coordinates. 
We used the numerical method, Eq. (44), to calculate the directional recoil rate dR/dcos9 
for SLI streams. The result is shown in Fig. 7 by the solid line. For comparison, we also 
show the result from the standard dark halo (dashed line), recomputed for the new reference 
direction. 

We see in Fig. 7 that dR/d cos 0|sli (solid line) peaks in the direction opposite to the 
case of dR/d cos 0| std (dashed line). This is because the average velocity of the SLI streams 
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FIG. 7: The directional recoil rate dR/dcos9 off a CS2 target for SLI streams (solid line) and 
the standard dark halo (dashed line) for a reference direction opposite to the direction of Galactic 
rotation, (A axis , p axis ) = (167?340, -59?574). 



points in a direction roughly opposite to that of standard dark halo's. Indeed, using Table 
I, the average (weighted) velocity Vsli of the SLI streams is 



V SLI = (57.534, -1.465, 85.994) km/s . 



(80) 



Comparison with Eq. (73) shows that V st d and Vsli form an angle of 170°. 



D. Anisotropic logarithmic-ellipsoidal models 

Many observations and numerical simulations show that galaxy halos are better approx- 
imated by triaxial models with anisotropic velocity distributions (e.g. Moore et al. [40], 
Helmi, White, and Springel [41], and Green [42]). Green [42] examined the effect of triaxial 
and anisotropic halo models on the exclusion limits from WIMP direct detection searches 
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and found that such models lead to non-negligible changes in the exclusion limits. Helmi, 
White, and Springel [41] found that the expected signal for the fastest moving DM particles 
in direct detection experiments is highly anisotropic. In this subsection we investigate the 
nuclei's directional recoil rate dR/dcos9 in specific cases of anisotropic logarithmic-ellipsoidal 
models. 

Evans, Carollo, and de Zeeuw [43] provided analytic solutions to the Jeans equations 
for the logarithmic-ellipsoidal model, which is the simplest triaxial generalization of the 
isothermal sphere, under the assumption of conical alignment of the velocity ellipsoid. Their 
anisotropic velocity distribution at the Sun's position can be approximated by an anisotropic 
Gaussian, which in the Galactic rest frame where the WIMPs are on average at rest is given 
by 

/(v) = , ,, 3/ * - - exp (S- ~ A- ~ A- V (81) 



(2 tt) 3 / 2 <j vX <t vY cr vZ V 2cr vX 2a vY 2a vZ/ 
Here X points toward the Galactic Center, Y toward the direction of Galactic rotation, 

and Z toward the North Galactic Pole. In this frame, the velocity dispersion matrix cr 2 is 

diagonal. 

In one of the Evans, Carollo, and de Zeeuw models, the Solar System is on the long 
(major) axis of the halo density ellipsoid, and 

v 2 

^ = (2 + 7)(p- 2 V 2 -l) ' (82) 
2 Vq(2p~ — 1) 

° vY = 2(p-* + q~ 2 - 1) ' (83) 



^o 2 (2g- 2 -i) 

vZ ' ' 2(p~ 2 + q~ 2 - 1) 



(84) 



In another of their models, the Solar System is on the intermediate (minor) axis of the halo 
density ellipsoid, and 

a '« x ~ (2 + 7 )(l + g- 2 -p" 2 )' (85j 



= 2(l + q - 2 -p- 2 ) ' (86) 

v 2 (2q~ 2 -p- 2 ) 
~ 2(l + q- 2 -p- 2 ^ ■ (87) 



Here p and q are constants used to describe the axis ratios of the density ellipsoid, and 7 
is a constant used to describe the velocity anisotropy. In the spherical limit (p = q = 1, 
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7 = 0), all the velocity dispersion components are equal, and the logarithmic-ellipsoidal 
model reduces to the isothermal sphere. 

An anisotropic Gaussian velocity distribution can be written in matrix form as 



/(v) = 



1 



■ cxp 



(2vr) 3 /2 (det <T2)l/2 

Its Radon transform has been found to be [23] 
f(w,w) 



--(v-V)V-V-V) 



(88) 



■ cxp 



[w - w • Vf 
2 w T crl w 



(89) 



(2ttw t cr^w) 1 ^ 

Since the WIMPs are assumed to be on average at rest in the Galactic frame, the same 
average WIMP velocity V as for the standard dark halo V st d (see Eq. (73)) applies to these 
anisotropic models. Since the velocity dispersion matrix cr^ i s diagonal in the Galactic frame, 
the principal axes of the velocity dispersion ellipsoid are aligned with the axes X, Y, and Z 
(see Eq. (81)). In matrix form, the exponent of Eq. (81) can be written as 



v x v Y v z 



-V o 


o\ 


l vx ) 


o -4- 

o v 

vY 





Vy 







\v z ) 



(90) 



and the velocity dispersion tensor is 



-1 = 



vX 







V 



a 2 vY 
al z) 



(91) 



To calculate w <r„w in Eq. (89), we write 



" T 2 ~ 

w a r w 



Wx Wy Wz 



(7 



vX 







a 2 vY 



V 



al z) 



wy 
ywzj 

Thus, the anisotropic recoil momentum distribution is given by 



a vX ™X + <Y ™Y 



°vZ W Z ■ 



(92) 



f(w,w) 



1 



■ exp 



[w - w • V] : 



V 2 *i°lx ™x + °Iy »y + °lz ™l) 
where only the components of a v and w are now in Galactic coordinates 



2 fe ™ 2 X + a lY ™ 2 Y + °lz ™ 2 z) J ' 

(93) 
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FIG. 8: The directional recoil rate dR/dcos 9 off a CS2 target as a function of cos 6 for logarithmic- 
ellipsoidal anisotropic models in which the Solar System is on the minor axis of the density ellipsoid 
(solid, dashed, dotted, and dashed-dotted lines, with parameters given in the legend). Also shown 
is the case of the isotropic standard dark halo (thick solid line). The reference direction n is 
opposite to the direction of the Galactic rotation, (A ax i s , /?axis) = (167?340, — 59?574). 

We fix the reference direction to be opposite to the direction of Galactic rotation, 
(^axis, fraxis) = (270°, 0°) in Galactic coordinates or (A axis , /3 axis ) = (167?340, -59?574) in the 
ecliptic coordinates. We consider eight anisotropic models (i.e. different values of p, q, and 
7) taken from Green [42]. Using the numerical method, Eq. (44), the results are shown in 
Fig. 8 in the case when the Solar System is on the minor axis of the density ellipsoid and 
in Fig. 9 in the case when the Solar System is on the major axis of the density ellipsoid. 
For comparison, in both figures we also show the case of the isotropic standard dark halo 
discussed in Sec. IV B (thick solid line). The values of p, q, and 7 in Figs. 8 and 9 are taken 
from Table I in [42]. 
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FIG. 9: Same as Fig. 8 but for logarithmic-ellipsoidal models in which the Solar System is on the 
major axis of the density ellipsoid. 

In Figs. 8 and 9, dR/dcos9 is an increasing function of cos6>. It is maximum in the 
forward direction (8 = 0°) because of the direction of the average WIMP velocity is similar 
to that of the standard halo's. We see from Figs. 8 and 9 that the behavior of the eight 
anisotropic models (solid, dashed, dotted, and dashed-dotted lines) resembles that of the 
isotropic standard dark halo (thick solid line). The differences between these models in 
Figs. 8 and 9 arise from the different velocity dispersions a v of these models resulting from 
different values of the parameter p, q, and 7 in Eqs. (82-87). This also means different ratios 
a v /V for each model. 

V. "FOLDED" DIRECTIONAL RECOIL RATE 



In this section, we will discuss the consequences of a lack of head-tail discrimination in 
WIMP direct detectors. We use the "folded" directional recoil rate dR/d | cos 9\ defined in 
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Eq. (9). We also discuss the exposures required to distinguish the expected WIMP signals 
from an isotropic noise, and present results for both an ideal zero-threshold detector and a 
detector with a finite energy threshold of 20 keV. We consider a CS 2 target (spin-independent 
and spin-dependent, which is zero) and a CF 4 target (assuming a spin-dependent cross 
section only). 

We use Eq. (9) and the numerical method, Eq. (44), for the standard dark halo (Sec. IV B), 
SLI streams (Sec. IV C), and logarithmic-ellipsoidal anisotropic models (Sec. IV D). The re- 
sults are shown in Figs. 10-13 for the CS 2 target nuclei, where for all cases the reference direc- 
tion is opposite to the direction of the Galactic rotation, (A ax i S , /3 ax i S ) = (167?340, — 59?574). 
In all of these figures, we show the rates for the standard dark halo and the SLI model. For 
the anisotropic models, the case in which the Solar System is assumed to be on the minor 
axis of the density ellipsoid is shown in Fig. 10 (for an ideal zero-threshold detector) and 11 
(for a 20-keV threshold detector). Figs. 12 and 13 show the analogous cases for the Solar 
System on the major axis of the density ellipsoid. 

The spin-dependent case for a 20-keV threshold detector is shown in Fig. 14 (when the 
Solar System is assumed to be on the minor axis of the density ellipsoid) and in Fig. 15 
(when the Solar System is assumed to be on the major axis of the density ellipsoid). The 
shape of the dR SD (> 20 keV)/d\cos9\ curves in Figs. 14 and 15 differs from the shape of the 
<iR(>20keV)/<i|cos#| curves in Figs. 11 and 13 respectively because of the different masses 
of sulfur (S) and fluorine (F). 

The dR/d\cos6\ curves that show small variation with \cos6\ will be harder to differentiate 
from an isotropic background. To make this statement quantitative, we compute the effective 
exposure £ required to distinguish a directional signal of WIMPs from a distribution uniform 
in | cos 6 1 that may be due to background events. For this purpose, we use the Kolmogorov- 
Smirnov test, which tests if data are drawn from a given distribution (in this case, the uniform 
distribution). Neglecting the background, we determine the required effective exposure £ as 
follows. We use a model dR/ d\ cos 9\ to Monte-Carlo generate a \co$6\ distribution of n events 
representing the outcome of a simulated experiment with zero background. We do this for 
10,000 simulated experiments. For each experiment, we compute the Kolmogorov-Smirnov 
statistic D n , which is the maximum vertical distance between the cumulative distributions of 
a uniform variate and of the simulated |cos#| values. We declare that the |cos#| distribution 
is non-uniform at the 90% level of significance (i.e. that the experiment under question 
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FIG. 10: The "folded" directional recoil rate dR/d\cos6\ off a CS2 target for the standard dark 
halo (thick solid line), SLI streams (thick dashed-dotted line), and four logarithmic-ellipsoidal 
anisotropic models in which the Solar System is on the minor axis of the density ellipsoid 
(solid, thick dotted, dashed, and dotted lines, with parameters given in the legend). In all 
cases, the reference direction is opposite to the direction of the Galactic rotation, (A ax is, Aixis) = 
(167?340,-59?574). 

detects a WIMP signal at the 90% significance level) when D n is greater than the upper 
10% quantile -D n l0 of the D n distribution evaluated under the null hypothesis of a uniform 
variate. To evaluate the probability distribution of D n under the null hypothesis, we use 
formulas in Refs. [44, 45] as implemented numerically in Refs. [33, 46]. We determine the 
fraction of the 10,000 simulated experiments with a positive detection at the 90% significance 
level. We finally increase the number n of events in each simulated experiment, until the 
fraction of experiments with a positive detection reaches 90%. This gives us the minimum 
number of events N e required to distinguish each of the halo models considered from a flat 
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FIG. 11: Same as Fig. 10 but using a detector with threshold energy equal to 20 keV. 

distribution at 90% significance level in 90% of the simulated experiments. Finally, the 
corresponding effective exposure £ is obtained by dividing N e by the integrated expected 
rate above threshold J c? | cos t9 1 [dR(> E thl )/d | cos 9\] (or the total rate J <i|cos [dR/d\cos9\] 
for an ideal zero-threshold detector). 

Table II shows the resulting N e and £ for a zero threshold detector and a detector with 
threshold energy equal to 20 keV, assuming a CS 2 target, a WIMP mass m = 60 GeV and 
a WIMP-proton cross section er p = 10~ 44 cm 2 . 

The number of events N e is independent of the value of the WIMP-proton cross section 
Op, because it depends only on the shape, and not the height, of the dR/d | cos 9\ distribution. 
The corresponding exposure however scales as the inverse of <t p . Regarding the dependence 
on the WIMP mass, in general N e and £ have a complicated dependence due to the relation 
between the threshold energy, the target nucleus mass, and the WIMP mass. 

From Table II, we notice the following. Even in the most favorable case of negligible 
background, the most optimistic case is the standard dark halo where, for a detector with a 
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FIG. 12: The "folded" directional recoil rate dR/d\cos 9\ off a CS2 target for the standard dark 
halo (thick solid line), SLI streams (thick dashed-dotted line), and four logarithmic-ellipsoidal 
anisotropic models in which the Solar System is on the major axis of the density ellipsoid 
(solid, thick dotted, dashed, and dotted lines, with parameters given in the legend). In all 
cases, the reference direction is opposite to the direction of the Galactic rotation, (A ax is, Aixis) = 
(167?340,-59?574). 

20-keV threshold, only 39 recoil events or an effective exposure of 113 kg-yr are needed in 
order to distinguish the standard dark halo from a flat distribution at 90% significance level 
in 90% of the simulated experiments. For a zero threshold detector, the required number 
of recoil events N e increases to 260 events (330 kg-yr of effective exposure). These numbers 
are almost the same as those obtained by Morgan, Green, and Spooner [26] for the same 
standard dark halo, although they used the average of \cos6\ while we use the full |cos#| 
distribution. 

In the case of SLI streams, the required number of events N e is 115 (295 kg-yr of effective 
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FIG. 13: Same as Fig. 12 but using a detector with threshold energy equal to 20 keV. 

exposure) for a detector with 20-keV threshold, and iV e =1606 (1596 kg-yr of effective expo- 
sure) for a zero threshold detector. Distinguishing an SLI streams signal from an isotropic 
background requires 3 to 5 times larger exposures than for a standard dark halo. 

The anisotropic models we considered are, with few exceptions, intermediate between the 
standard dark halo and the SLI streams cases. The hardest cases to distinguish from a flat 
distribution are the two anisotropic models with parameters p = 0.72, q = 0.7, 7 = 4.02 and 
p = 0.72, q = 0.7,7 — 2.01. They require effective exposures of thousands to ten thousands 
of kg-yr. This is due to the peculiar behavior of their <iR(>20keV)/<i|cos#| distributions. 
As seen in Fig. 11, these distributions (dashed and dotted lines) exhibit a sudden drop in 
dR(> 20 keV) /d\ cos 9 | in the folded forward+backward direction | cos t9| = 1. Notice that 
these two models are the most anisotropic models among the eight anisotropic logarithmic- 
ellipsoidal models considered in this study (7 = 4.02 and 7 = 2.01, respectively). This might 
be the reason for such behavior. 

To illustrate a spin-dependent case, we assume a CF 4 target, a proton-odd approximation 
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FIG. 14: The spin-dependent "folded" directional recoil rate dR SD /d\cos 6\ off a CF4 target using 
a detector with threshold energy equal to 20 keV, for the standard dark halo (thick solid line) , SLI 
streams (thick dashed-dotted line), and four logarithmic-ellipsoidal anisotropic models in which 
the Solar System is on the minor axis of the density ellipsoid (solid, thick dotted, dashed, and 
dotted lines, with parameters given in the legend). In all cases, the reference direction is opposite 
to the direction of the Galactic rotation, (A axis , /3 axis ) = (167?340, -59?574). 

for the F spin, a spin-dependent cross section off protons a^ D = lCT 44 cm 2 , and a vanishing 
spin-independent cross section. We examine both a zero threshold detector and a detector 
with threshold energy equal to 20 keV. Because the spin- dependent rates do not increase 
as A 2 , the effective exposures required for spin-dependent interactions are longer than for 
the spin-independent case (at the same WIMP-proton cross section). For CF 4 we find them 
about 1,000 times longer than for CS 2 . For example, in the case of the standard dark halo 
we find that £ is about 4.8xl0 5 kg-yr for a zero threshold detector and it is about 1.6xl0 5 
kg-yr for a detector with threshold energy equal to 20 keV. The effective exposures £ for 
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FIG. 15: Same as Fig. 14 but for logarithmic-ellipsoidal models in which the Solar System is on 
the major axis of the density ellipsoid. 



the other halo models follow the same pattern and they are in the order of magnitude of 10 6 
kg-yr. Such effective exposures are impractical. 

To summarize, the "folded" directional recoil rate dR/d\cos$\ can be helpful in recogniz- 
ing the cases of the standard dark halo, SLI streams, and some not-too-anisotropic models. 
However, if the detector threshold is too low, or the degree of anisotropy too high, it may 
be difficult to recognize SLI streams and some other anisotropic models. 



VI. SUMMARY AND CONCLUSIONS 



In this paper, we studied the directional recoil rate dR/dcos8 of recoiling target nuclei 
struck by WIMPs in terms of the angle 9 between the nucleus recoil direction w and a chosen 
reference direction n in the sky. We used the ecliptic coordinate system and imagined a CS2 
detector similar to DRIFT but with 3D read-out capabilities. 
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TABLE II: Number of recoil events N e required to distinguish each of the halo models considered 
from a flat distribution at 90% significance level in 90% of the simulated experiments using a 



Kolmogorov-Smirnov statistic, for a CS2 zero-threshold detector and a CS2 


detector with threshold 


energy equal to 20 keV. Here a p = 10 44 cm 2 . Also shown 


are the effective exposures £ required 


for each case. 






zero threshold energy 

Model 




threshold energy=20 keV 


N e £ (kg-yr) 




iV e £ (kg-yr) 


standard dark halo 260 330 




39 113 


SLI streams 1606 1596 




115 295 


Anisotropic logarithmic-ellipsoidal models (Solar System is on the minor 


axis of the density ellipsoid): 


p = 0.9, q = 0.8, 7 = 0.07 503 585 




98 275 


p = 0.9, q = 0.8,7 = -0-62 876 969 




196 482 


p = 0.72, g = 0.7, 7 = 4.02 1005 1149 




1199 3548 


p = 0.72, q = 0.7,7 = 2.01 2088 2276 




7500 19116 


Anisotropic logarithmic-ellipsoidal models (Solar System is on the major axis of the density ellipsoid): 


p = 0.9, q = 0.8,7 = -1 365 427 




60 168 


p = 0.9, q = 0.8, 7 = -1.33 600 668 




115 283 


p = 0.72, q = 0.7, 7 = -1.39 301 354 




46 130 


p = 0.72, q = 0.7, 7 = -1.6 497 555 




90 221 


The directional recoil rate dR/dcos$ was computed and compared for different halo mod- 


els that represent several WIMP velocity distributions: 


streams of WIMPs (Sec. IV A), the 



standard dark halo (Sec. IV B), Sikivie's late-infall halo model (SLI streams) (Sec. IV C), 
and anisotropic logarithmic-ellipsoidal models (Sec. IV D). We repeated our analysis for a 
"folded" directional recoil rate dR/d\cos6\ that incorporates the inability of some detectors 
to distinguish the beginning of a recoil track from its end (lack of head-tail discrimination). 
For all of the halo models considered, we compared dR/d\cos8\ to an isotropic background, 
to examine the possibility of discriminating a WIMP signal from background noise. 

We computed dR/dcos6 both numerically and analytically (Sec. III). The numerical 
method (Sec. Ill A) uses a fifth-order Cash-Karp Runge-Kutta method and can be applied 
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to general (Gaussian and non-Gaussian) velocity distributions and any reference direction. 
The analytical method (Sec. Ill B) works only for Gaussian distributions with the reference 
direction n aligned with the average WIMP velocity V. The analytic formula was used to 
cross check the numerical calculation. Comparison of the numerical and analytic calculations 
gave the same results. In both the numerical and analytical methods, the recoil momentum 
function f(w,w) used in the calculation of dR/dcos9 was taken as the Radon transform of 
the velocity distribution function /(v). 

We applied the numerical method to the aforementioned WIMP halo models and pre- 
sented the results by showing the directional recoil rate dR/dcos9 as a function of the angle 
9 (see Sec. IV). For generic streams of WIMPs, we showed how varying the ratio cr v /V of 
the velocity dispersion a v to the magnitude of the average WIMP velocity V affect the direc- 
tional recoil rate dR/dcos9 (see Fig. 4). We also showed the effect of varying the reference 
direction n or equivalently the stream velocity V (see Fig. 5). 

Comparisons between the case of the SLI streams and the case of the standard dark halo 
showed that SLI streams produce a directional recoil rate that peaks in the opposite direction 
to the standard halo one (see Fig. 7). The case of streams with anisotropic logarithmic- 
ellipsoidal models resembles that of the standard dark halo, with small differences between 
the anisotropic models due to different values of their axial ratios p and q and anisotropy 
parameter 7 (see Figs. 8 and 9). 

We allowed for the difficulty of head-tail discrimination in WIMP direct detection exper- 
iments in Sec. V. There we introduced a "folded" directional recoil rate dR/d\cos9\ suitable 
for direct comparison with experiments lacking head-tail discrimination. For both a zero 
threshold detector and a detector with 20-keV threshold energy, we calculated the number of 
recoil events N e and the effective exposure £ required to distinguish each of the halo models 
considered from a flat distribution (see Table II). We found that in distinguishing a signal 
from an isotropic background noise, the "folded" directional recoil rate dR/d\cos9\ may be 
effective for the standard dark halo and some of the anisotropic logarithmic-ellipsoidal mod- 
els; it may be less effective for the SLI streams and other anisotropic models (see Figs. 10- 
13). In most cases, for m = 60 GeV and a v = 10~ 44 cm 2 , exposures from few dozens to 
few hundreds of kg-yr of CS2 would be needed to utilize the "folded" directional recoil rate 
dR/d\cos9\ for the purpose of discriminating a directional WIMP signal from an isotropic 
background noise. 
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